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Nucleon electromagnetic form factors are fundamental quantities related to the charge 
and magnetization distributions inside the nucleon. Understanding the nucleon electro- 
magnetic structure in terms of the underlying quark and gluon degrees of freedom of 
quantum chromodynamics is a challenging and urgent task. The nucleon electromag- 
netic form factors have been studied in the past extensively from unpolarized electron 
scattering experiments. With the development in polarized beam, recoil polarimetry, 
and polarized target technologies, polarization experiments have provided more precise 
data on these quantities. At the same time, significant theoretical progress in areas rang- 
ing from effective field theories to lattice QCD calculations, has been made in describing 
these data. In this article, I review recent experimental and theoretical progress on this 
subject. I will also provide future outlook on this topic. 

1. Introduction 

Nucleons (protons and neutrons) arc fundamental building blocks of matter. 
Nevertheless, they are not fundamental particles. Otto Stern received the Nobel 
Prize in Physics in 1943, for his experimental methods of studying the magnetic 
properties of nuclei, in particular for measuring the magnetic moment of the proton 
itself. The anomalous magnetic moments of the proton and neutron reveal that 
nucleons are not point-like Dirac particles, but particles with underlying structure. 

The electromagnetic form factors of the nucleon have been a longstanding sub- 
ject of interest in nuclear and particle physics. They are fundamental quantities 
describing the distribution of charge and magnetization within nucleons. Probing 
the nucleon electromagnetic structure has been an ongoing experimental endeavor 
since the discovery of the anomalous magnetic moment of the proton. 

Electron scattering has been proven to be a very useful tool in probing the struc- 
ture of nucleon and nuclei. The leptonic part of the vertex is well understood from 
Quantum Electrodynamics (QED), thus it is a clean probe of hadronic structure. 
Furthermore, the electromagnetic coupling constant is relatively weak, so higher 
order diagrams are suppressed compared to the lowest order one-photon-exchange 
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diagram. The study of the electromagnetic structure of the proton from electron- 
proton elastic scattering was pioneered by Hofstadter and colleagues at the Stanford 
Linear Accelerator Center (SLAC) in the 1950s, for which Hoftstadter received the 
Nobel Prize in Physics in 1961. The deep-inelastic scattering experiments of elec- 
trons from protons carried out by Friedman, Kendall and Taylor which established 
the underlying quark structure of the proton at SLAC, led to the Nobel Prize in 
Physics in 1990. 

Quantum Chromodynamics (QCD) is the theory of strong interaction in terms 
of quark and gluon degrees of freedom. While QCD has been extremely well tested 
in the high energy regime, where perturbative QCD is applicable, understanding 
confinement and hadron structure in the non-perturbative region of QCD remains 
challenging. Knowledge of the internal structure of protons and neutrons in terms 
of quark and gluon degrees of freedom is not only essential for testing QCD in the 
confinement regime, but it also provides a basis for understanding more complex, 
strongly interacting matter at the level of quarks and gluons. 

The rest of the paper is organized as follows: Section II gives an introduction on 
the nucleon electromagnetic form factors; Section III contains a discussion on form 
factor data from unpolarized electron scattering experiments; Section IV presents 
recent data on the nucleon electromagnetic form factors from double-polarization 
experiments; Section V reviews theoretical progress on the nucleon electromagnetic 
form factors; and Section VI provides future outlook on this subject. 

2. The Nucleon Electromagnetic Form Factors 

From QED, the lowest-order amplitude (Fig. 1) for electron-nucleon elastic scat- 
tering is given by 

Tf, = -ij 3,{^)J^d^x, (1) 
where q = p' — p and the electron transition current is 

f = -e«(fc')7''u(A;)e*('='-'=)-^ (2) 

The nucleon is an extended spin-i object, thus the nucleon transition current is 
more complicated than that of the electron. Based on the requirements of covariance 
under the improper Lorentz group, current conservation and parity conservation, 
the nucleon transition current is written as 

?i(p)e^(f'-f)■^ (3) 

where Fi and F2 are two independent form factors, also called the Dirac and 
Pauli form factors, respectively, k is the anomalous magnetic moment, and M is 
the nucleon mass. 

Ernst, Sachs, and Wall ^ defined the following, the so-called Sachs' form factors: 
Ge{q^) (electric form factor) and GM{<f) (magnetic form factor), which are written 



J" = eu{p') F,iq^)r + ^F^il'y^'^l-^ 
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as linear combinations of Fi and F2: 




Fig. 1. The one-photon-exchange diagram for electron-nucleon elastic scattering. 

The Sachs' form factors, Ge and Gm, have more intuitive physical interpreta- 
tions than Fi and F2. In electron scattering, the inverse of the four-momentum 
transfer squared, {Q^ = —q^ > 0) is commonly used. In the limit of = 0, 

G^(g2 = 0) = l,G^(Q2 = o) = 

and 

= 0) = Mp, Gm(Q' = 0) = Mn, 

corresponding to the charge and the magnetic moment of the nucleon, respectively. 

Sachs ^ demonstrated that Ge and Gm can be interpreted as Fourier transforms 
of spatial distributions of charge and magnetization of the nucleon in the so-called 
Breit frame. For elastic electron-nucleon scattering, the Breit frame coincides with 
the center-of-mass frame of the electron-nucleon system. In this reference frame, 
the incoming c>lcK;tron has a momentum of | and the nucleon initial momentum 
is — I ; the scattered electron has a momentum of — | and the recoil proton has a 
momentum of |. Thus, the Breit frame is a special Lorentz frame in which = — ^, 
i.e., no energy transfer is involved in this particular reference frame. Thus, for 
each q^ value, there is a Breit frame in which the form factors are represented as 
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GE,M{'f) = Ge,m{q^), where GE,M{q^) is determined in the lab frame. One can 
therefore perform a three-dimensional Fourier transformation once the form factor 
information is available, 

.«=/ge--^G.(,-). (6) 

Although relativity obscures the physical interpretation of such a quantity, it is 
analogous to the "classical" charge density distribution. 

In the low region, below the two-pion production threshold (i.e., q^ < to, 
where to = (2m,r)^ < 2 fm~-^) the energy transfer in the scattering process is 
negligible and nucleon electric and magnetic form factors can be taken as the Fourier 
transforms of the charge and magnetization radial distributions pchg{r) and Pmag('') 
in the rest frame of the nucleon. 

In the non-relativistic limit, the root-mean-square (rms) charge radius of the 
nucleon is related to the nucleon electric form factor as: 

<r^> dGEm , 

with the boundary condition Ge{Q^ = 0) = 1 for the proton and for the neutron. 
The corresponding rms magnetic radius is 

<r^> 1 dGMm. 

where /i^v is the nucleon magnetic moment, /Xp =2.79 and /U„ =-1.91, in the units 
of the nucleon magneton. 

Recent results from lattice QCD calculations suggest that the nucleon rms 
charge radius can be calculated from first principles with an uncertainty of only a few 

percent, and this field is rapidly evolving due to both improvements in compiiter 
architecture and new algorithms. Thus, precise information on this fundamental 
quantity is essential in terms of testing the QCD prediction from the lattice. 

Accurate information about the proton charge radius is also essential in con- 
ducting high-precision tests of QED from hydrogen Lamb shift measurements. The 
standard Lamb shift measurement probes the 1057 MHz fine structure transition 
between the 2Si/'2 and 2Pi/2 states in hydrogen. The hydrogen Lamb shift can be 
calculated to high precision from QED using higher order corrections. The proton 
rms charge radius is an important input in calculating the hadronic contribution to 
the hydrogen Lamb shift. 

The two most precise and widely cited determinations of the proton charge 
radius from electron scattering experiments in the literature give Vp = 0.805(11) 
fm ^ and Vp = 0.862(12) fm ^, respectively, differing from each other by more than 
7%. While the recent precision hydrogen Lamb shift measurements 7,8,9, io,ii g^j.^ 
better agreement with the QED predictions using the smaller value of the proton 
charge radius without the two-loop binding effects, they are consistent with the 
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larger value of the proton charge radius when two-loop binding effects are included in 
the QED calculations. Mostly recently, a QED calculation with three-loop binding 
effects has been carried out Before accurate comparisons between theory and 
experiment can be made, either in QCD or QED, a new, high precision experiment 
on the proton charge radius is urgently needed. 

The mean squared charge radius of the neutron is < >= — 0.113±0.003±0.004 
fm^, obtained most recently from an experiment carried out at the Oak Ridge 
National Laboratory in which the neutron transmission through a thorogenic 
liquid ^°^Pb has been measured in the neutron energy range between 0.1 and 360 



Among the four nucleon electromagnetic form factors, the neutron electric form 
factor, is of particular interest, but it is also least known due to the lack of free 
neutron targets and the smallness of in general. Its interpretation within many 
models has been obscured by relativistic effects. Recently, Isgur demonstrated 
to leading order in the relativistic expansion of a constituent quark model that 
the naive interpretation of G^ as arising from the neutron's rest frame charge 
distribution holds due to a cancellation of the Foldy term against a contribution to 
the Dirac form factor, Fi . Enormous experimental progress has been made on the 
nucleon electromagnetic form factors in the last decade or so. In particular, new 
data with significantly improved precision from double-polarization experiments are 
available due to the recent advances in polarized target, beam and recoil polarimeter 
technologies. More precise data from the next generation of double-polarization 
experiments are anticipated in the next several years, and these data will allow more 
stringent tests of theoretical descriptions of the nucleon electromagnetic structure. 
In the next two sections, the experimental progress on the nucleon electromagnetic 
form factors will be reviewed, in particular the most recent measurements from 
double-polarization experiments. 

3. Unpolarized Electron Scattering and Nucleon Electromagnetic Form 
Factors 

3.1. Elastic Electron- Proton Scattering 

In terms of G^ and Gm, the differential cross-section for electron- nucleon scat- 
tering can be written in the one-photon-exchange picture as: 



where E' and 6 are the scattered electron energy and angle, respectively; a is the 
fine structure constant, and r = -^jj- 

The proton electric (G^) and magnetic (G^) form factors have been studied 

extensively in the past from unpolarized electron-proton (ep) elastic scattering using 
the Rosenbluth separation technique Eqn. 9 can be re- written as: 



eV. 
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da 



<J"E H <J"n, 



1 + T 



Cm 



1 



(10) 



where e = (1 + 2(1 + r) tan 



2 e^ 



is the virtual photon longitudinal polarization, 



and ctm is the Mott cross section describing the scattering from a pointlike target: 

a2i;'cos2| 

CM = , ■ 4 ■ (11) 
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Fig. 2. World data published since 1970 on the proton electric form factor obtained by the Rosen- 
bluth method from the unpolarized cross-section measurements. The data are plotted as q^, 
where Go is the standard dipole parameterization. 



In the Roscnbluth method, the separation of and is achieved by 

measuring the cross section at a given value by varying the incident electron 
beam energy and the electron scattering angle. The measured differential cross 
section is then plotted as a function of tan^ |, and one can extract information 
on G^'^ and G^^ from the slope and the intercept of the plotted curve. While 
the G^^ term dominates the cross section in the low region, the G^^^ term 
dominates at large values. Thus, the extraction of G^ at low Q^, and G^ 
at large values becomes difficult using the Roscnbluth technique. Fig. 2 and 
Fig. 3 show the world data since 1970 on the proton electric and magnetic form 
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factors, as a function of using the Rosenbluth separation technique. The data 

are shown as ^j^, and ^ , respectively, where /ip is the proton magnetic moment, 
and Gd = (1 + Q^/O.Tl)"^ is the standard dipole parameterization, and is in 
(GeV/c)2. 

3.2. Electron- Deuteron Scattering 

Because of the lack of free neutron targets, little is known about the neutron 
electromagnetic structure. The best known quantity is the neutron rms charge ra- 
dius which is obtained from thermal neutron scattering experiments. The neutron 
electromagnetic form factors are known with much less precision than the proton 
electric and magnetic form factors. They have been deduced in the past from elastic 
or quasielastic electron-deuteron scattering. This procedure involves considerable 
model dependence. Another complication arises from the fact that the net charge 
of the neutron is zero. As such the neutron electric form factor is much smaller 
than its magnetic form factor G^. Therefore, the magnetic part of the contribu- 
tion dominates the cross section, which makes it very difficult to extract from 
unpolarized cross section measurements using deuterium targets. 




(GeV/c) 



Fig. 3. World data published since 1970 on the proton magnetic form factor obtained using the 

qP 

Rosenbluth method from unpolarized cross-section measurements. The data are plotted as 



in the unit of the standard dipole parameterization. 



3.2.1. The Elastic Electron- Deuteron Scattering 

The cross section for the unpolarized elastic electron-deuteron scattering in the 
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one-photon-exchange approximation is described by the Rosenbluth formula, 



da 



(12) 



The two structure functions: A{Q'^) and B{Q'^) can be separated by the Rosenbluth 
separation technique discussed previously. The deuteron is a spin-1 nucleus and the 
characterization of its charge and magnetization distribution requires three form 
factors: Fo(Q^), -Fq(Q^), -FmCQ^), the charge monopole, quadrupole and the mag- 
netic dipolc form factor, respectively. The structure functions Ai^Q"^) and B{Q^) 
can be expressed in terms of these three form factors as 



(13) 



-t{1 + t)FUQ% 



(14) 



where r 



4Mj 



Thus, it is not possible to separate all three form factors of 



the deuteron from the unpolarized elastic electron-deuteron cross section measure- 
ment alone. An additional measurement which involves polarization is necessary 
to separate the charge monopole and the charge quadrupole form factors. Such a 
polarization experiment can be either a deuteron tensor polarization measurement 
employing a polarimeter by using an unpolarized electron beam and an unpolarized 
deuteron target or an analyzing power measurement by using a polarized deuteron 
target and an unpolarized electron beam. 

The tensor moment t2o is particularly interesting due to its sensitivity to Fc- It 
is defined as 



1 



< - 2 > = — 



1 

V2 



N++N_- 2No 



N+ + N_+ No 



(15) 



where N^, and Nq are the occupation numbers for the deuteron magnetic sub- 
state Ms = 1, 0, —1, respectively. The deuteron tensor moment t2o can be expressed 
in terms of the deuteron form factors as: 



t20 



-1 



V2I0 



-tFcFq + -r^F^ + -r(l + 2(1 + r) tan^ -)Fj 



M 



(16) 



where Iq = A{Q^) + B(Q^) tan"^ |. Therefore, by combining the structure functions 
A{Q'^), B{Q'^) from the unpolarized cross section measurement, and the deuteron 
tensor moment ^20 measurement, one can separate all three deuteron form factors. 

The structure function AiQ"^) provides one of the few methods to infer the neu- 
tron electric form factor, especially in the low region (less than 1.0 (GeV/c)^) 
where theoretical descriptions of A(Q^) including relativity, meson-exchange cur- 
rents (MEG), etc. are under better control compared to higher region. The most 
systematic information on at low Q^, prior to any polarization experiment, is 
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from the A{Q^) structure function determined from the elastic electron-deuteron 
scattering experiment by Platchkov et al. However, the extraction procedure 
is quite compHcated. First, the subtraction of F^{Q^) from A{Q^) using data on 
B{Q^) is performed to obtain the corrected A{Q'^) which contains contributions 
from Fc{Q^), and Fq{Q'^) only. Second, the rclativistic and MEC corrections arc 
apphed to the corrected A{Q'^) to obtain the corresponding A{Q'^) in the impulse 
picture. Next, the deuteron structure is removed to obtain the nucleon iso-scalar 
charge form factor. Finally, the proton electric form factor is subtracted from the 
nucleon iso-scalar charge form factor and is obtained. The extracted values 
are extremely sensitive to the deuteron structure. Fig. 4 shows the G^ values ex- 
tracted with the Paris nucleon-nucleon potential together with a fit of the data 
(dash-dotted curve). Results from fitting the G^ data extracted with the Nijmegen 
potential the Argonne V14 (AV14) and the Reid-Soft Core (RSC) NN 
potentials are shown as solid, dashed and dotted curves, respectively. The large 
spread represents the uncertainty of G^ due to the deuteron structure, and the 
absolute scale of contains a systematic uncertainty of about 50% from such an 
extraction. 



nijmegen 
Av14 
peins 
rsc 




O 0.2 0.4 0.6 

Fig. 4. The electric form factor of the neutron, G^, as a function of four-momentum transfer 
squared from the unpolarized electron-deuteron elastic scattering measurement by Platchkov et 
al. 

Recently, Schivilla and Sick extracted G^ from an analysis of the deuteron 
quadrupole form factor Fq{Q^) data. Such an approach is different from the anal- 
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ysis employing the structure function A{Q'^) discussed above, in which both the 
charge monopole and quadrupole deuteron form factors contribute. State-of-the- 
art deuteron calculations based on a variety of different model interactions and 
currents show that the Fq form factor is relatively insensitive to the uncertain two- 
body operators of shorter range because the long-range one-pion exchange operator 
dominates the two-body contribution to Fq . Also, data on the Fq form factor has 
improved tremendously with the measurements of the polarization observable t2o in 
electron-dcniteron elastic scattering. As such, the neutron electric form factor has 
been extracted up to a value above 40 fm^^ (1.55 (GeV/c)^) from Fq{Q'^) with- 
out undue systematic uncertainties from theory. These extracted values of will 
be presented later together with data from double-polarization experiments. 

3.2.2. Quasielastic Electron- Deuteron Scattering 

Quasielastic electron-deuteron scattering, in which the kinematics of the elec- 
tron scattering from the nucleon inside the deuteron is selected, is the other process 
involving the deuteron which has been used extensively in probing the electromag- 
netic structure of the neutron. It includes both inclusive measurements, in which 
only scattered electrons are detected at the quasielastic kinematics, and coinci- 
dence measurements where both the scattered electron and the knockout neutron 
are measured. 

The missing mass squared for quasielastic e — d scattering is defined as, = 
M2 + 2M{E - E') - where M is the nucleon mass; W'^ = = 0.88 GeV^ at 
the quasielastic peak. The measured quasielastic e — d cross-section per nucleon, 
<t{E, £", 6) converted to the reduced cross section is written as: 

aR = e{l + T' f^^'^''^^ =RT + eRL, (17) 

where <tm is the Mott cross section defined earlier, and Rt and R^ are the transverse 
and longitudinal nuclear response functions, respectively. In the plane wave impulse 
approximation (PWIA) the quasielastic Rt response function is proportional to 
Gm^ + G^^ and the Rl response function is proportional to G^"^ + . Thus, 
the extraction of the neutron electromagnetic form factor requires the separation of 
the Rl and the Rt response functions using the Rosenbluth technique, followed by 
the subtraction of the proton contribution in PWIA. 

Until recently, most data on G"^ had been deduced from quasi-elastic electron- 
deuteron scattering. For inclusive measurements 23,24,25,26,27^ ^j^g procedure requires 
the separation of the longitudinal and transverse cross sections and the subsequent 
subtraction of a large proton contribution. Thus, it suffers from large theoretical 
uncertainties due in part to the deuteron model employed and in part to corrections 
for final-state interactions (FSI), MEG effects, and relativistic corrections. The pro- 
ton subtraction can be avoided by measuring the neutron in coincidence d{e,e'n), 
but the difficulty is associated with the absolute neutron detection efficiency cali- 
bration. This method was used by Stein et al. Bartel et al. and most recently 
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at MIT-Bates Laboratory by Markowitz et al. . In addition, the anti-coincident 
d(e, e'p) measurement, where the absence of the proton detection is required, was 
also carried out in the past ^1,32 

The sensitivity to nuclear structure can be greatly reduced by taking the cross- 
section ratio of d{e,e'n) to d{e,e'p) at quasi-clastic kinematics. Several recent ex- 
periments 33,34,35,36 ji^ye employed this technique to extract with uncertainties 
of <2% 35,36 foj. Q2.^a,lues from 0.1 to 0.8 (GeV/c)^ While such level of precision 
is excellent, there is considerable disagreement among the results 30,33,34,35,36 ^^^^ 
cerning the absolute value of G^f- coincidence experiments require an absolute 
calibration of the neutron detection efficiency, which is a challenging task. Most 
recently, the measurement using the ratio technique has been extended to a 
value of above 4.0 (GeV/c)^ at Jefferson Lab While discrepancies among the 
deuterium experiments described above may be understood additional data on 
G%i , preferably obtained using a complementary method, are highly desirable. The 
inclusive quasi-elastic reaction ^lie{e, e') provides just such an alternative approach. 
In comparison to deuterium experiments, this technique employs a different target 
and relies on polarization degrees of freedom, and is therefore subject to completely 
different systematics. We will discuss this type of double-polarization measure- 
ments in the next section. The world data on GjJ^ from unpolarized ed quasielastic 
scattering experiments are shown in Fig. 5. 
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Fig. 5. The published neutron magnetic form factor data from unpolarized deuterium experiments 
[23-36]. The data are plotted as G'^fiXn, in units of the standard dipole form factor Go, as a 
function of . 
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Polarization degrees of freedom in electron scattering have proven to be very 
useful in extracting information about small amplitude processes by isolating terms 
sensitive to the interference of the small amplitude with a much larger amplitude. 
The ability to selectively isolate certain combinations of amplitudes that either 
polarized electron beams or polarized targets (recoil polarimcters) have historically 
provided when used in isolation, is significantly enhanced by using them together. 
Thus, polarized electron scattering from polarized targets, and polarization transfer 
measurements using recoil polarimeters can provide important information on the 
nucleon spin structure and the electromagnetic structure of the nucleon. 

The development of polarized beams, polarized targets and recoil polarimeters 
has allowed more complete studies of electromagnetic structure than has been pos- 
sible with unpolarized reactions. In quasielastic scattering, the spin degrees of 
freedom introduce new response functions into the differential cross section, which 
provides additional information on nuclear structure Such studies are especially 
important for probing the neutron electromagnetic structure because of the lack of 
free neutron targets in nature. Section 4.1 contains a brief overview of experimental 
techniques employed in double polarization experiments. 

4.1. Experimental Techniques in Double- Polarization Experiments 

4.1.1. Polarized electron beam 

Polarized electrons can be produced by various techniques Photoemission 
from GaAs has several advantages: high peak current, helicity reversal by optical 
means, and the feasibility of building a GaAs source with small transverse phase 
space or emittance, an important requirement for some electron accelerator ma- 
chines. 

The GaAs polarized electron source works through photo-emission of electrons 
that have been polarized through optical pumping. The electron beam is generated 
by illuminating a GaAs cathode, which is placed in ultra-high vacuum, with high- 
intensity circularly polarized laser light. Strained GaAs cathodes are frequently used 
nowadays and have advantage over conventional bulk GaAs cathodes. The strain 
creates a gap in the different sub-levels of the P3/2 electrons in the valence bands 
of the GaAs Left circularly polarized light with the right frequency, incident 
on the GaAs, only drives transition between the P3/2 m=3/2 state in the valence 
band and the 5*1/2 1/2 level of the conduction band. From there the polarized 
electrons diffuse to the surface and escape into the surrounding vacuum. Typically, 
the surface of the GaAs cathode is treated with cesium to create a negative electron 
affinity. The polarization of source electrons is consequently as high as 80%. 

The electrons from the polarized source are longitudinally polarized, and are 
then accelerated in the linac. The electron beam polarization can be determined 
either by using a M0ller polarimeter which measures the asymmetry from polarized 
electrons scattering off polarized atomic electrons in a magnetized iron foil, or by 
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using a Compton polarimeter in which the asymmetry from circularly polarized 
laser light backscattering from the polarized electron beam is measured. In some of 
the high precision asymmetry measurements, both polarimeters might be employed 
to minimize systematic uncertainties in the determination of the electron beam 
polarization. 

4.1.2. Polarized targets 

Polarized targets for studies of nucleon structure can be categorized as polarized 
proton targets and polarized effective neutron targets because of the lack of free, 
stable neutrons in nature. For the purpose of a brief overview of polarized targets, 
the discussion below is organized according to gaseous and solid targets. 

There are two general types of polarized gaseous targets: external and internal. 
External targets arc scaled targets with windows. They are miich more d(^nsc than 
internal targets in order to achieve a practical luminosity with relatively low inci- 
dent electron beam currents. The windowless internal targets placed in an electron 
storage ring with very high electron beam currents are open flowing systems, thus 
with much thinner average target thicknesses compared to external targets. 

The internal gas targets are typically polarized by either optical pumping or 
by the atomic beam source method. Optical pumping is widely used to polarize 
a sample of atoms by transferring angular momentum from a pump light beam, 
typically a laser beam, to the sample atoms. The atomic beam source (ABS) method 
is based on the well-known principle of Stern-Gerlach separation and radio frequency 
(RF) transitions. One can employ either the atomic beam source method which is 
a well-established technique, or the optical pumping method for polarized hydrogen 
(proton) and deuterium (effective neutron) internal targets. The typical atomic flux 
an ABS can feed a storage target cell is ~ 6 x 10^^ atoms/second and the proton 
and deuteron nuclear vector polarization can be as high as 90%. 

An internal polarized H/D gas target is also feasible with spin-exchange optical 
pumping method. Such a target works in the following way: polarized laser light 
optically pumps alkali atoms. Hydrogen (deuterium) atoms are polarized through 
spin-exchange collisions with the alkali atoms. At large hydrogen (deuterium) den- 
sities, frequent H-H (D-D) collisions enhance the probability of hyperfine interac- 
tions, and the system approaches spin temperature equilibrium The nucleus is 
polarized in spin temperature equilibrium, in which hydrogen nuclear and atomic 
polarizations are the same and, in the case of deuterium nuclear vector polarization 
exceeds the atomic polarization. Spin-exchange optical pumping of alkali-metal 
atoms with hydrogen (deuterium) is an efiicient way to polarize hydrogen (deu- 
terium) atoms and atomic fluxes as high as 10"'^* atoms/second are typical with this 
type of source, though with much lower nuclear polarization than that of ABS, a 
50% nuclear vector polarization for hydrogen and 60% for deuterium can be antic- 
ipated. 

Polarized '^He is useful for studying the neutron structure because its groimd 
state is dominated by a spatially symmetric S wave in which the proton spins cancel 
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and the spin of the ^He nucleus is predominantly carried by the unpaired neutron 
alone. Internal polarized '^He gas targets can be produced in principle either by 
using optical pumping methods or atomic beam source method. However, only 
metastability-exchange optical pumping has been used up to now in producing in- 
ternal targets of practical thicknesses for scattering experiments. In the case of '^He, 
direct optical pumping between its ground state and the first excited state is not pos- 
sible because of the energy difference involved. The metastability-exchange optical 
pumping involves optical pumping of 2^5'i metastable state atoms, then transfer- 
ring the polarization to ^He ground state atoms through metastability-exchange 
collisions. 

The only external gaseous targets that exist for scattering experiments are po- 
larized ^Hc targets based on both the metastability-exchange optical pumping tech- 
nique and the spin-exchange optical pumping technique in which rubidium vapor 
is typically used as the spin-exchange medium. While spin-exchange optical pump- 
ing is capable of producing high density targets, metastability-exchange optical 
pumping combined with mechanical compression produce polarized ^He targets of 
comparable densities. External polarized ^He targets based on these techniques 
have been used successfully at electron accelerator laboratories around the world. 
For more details about these gaseous polarized targets, we refer the reader to a 
review article by Chupp, Holt and Milner 

Polarized solid targets have been used widely in lepton scattering experiments 
probing nucleon structure, particularly in deep-in-elastic lepton scattering exper- 
iments studying nucleon spin structure. The SMC experiment at CERN used 
a low intensity muon beam incident on a very thick target of butanol or deuter- 
ated butanol, which was cooled by a powerful dilution refrigerator. Polarization is 
achieved by the technique of dynamic nuclear polarization, which works in the fol- 
lowing way. A hydrogeneous compoimd is doped with a small concentration (10^^) 
of unsaturated electron spins. In a large magnetic field and low temperatures these 
electrons are polarized from the Boltzmann distribution. At a temperature of 1 K 
and a magnetic field of 5 Tesla, the electron polarization is 99.8%. The relaxation 
time of the electrons in the higher-energy state is short (~ 10~^ second) due to 
the interaction of the electrons with the lattice. The electrons therefore flip back 
and are available for further spin-flip transitions with other protons (deuterons). 
The relaxation time of the nuclei, on the other hand, is long lO'' second). This 
allows the nuclei to accumulate preferentially in a state selected by the frequency 
of the microwaves, thus leading to a high polarization of the sample. Experiments 
at SLAC (E143, E155) ^^'^^ and Jefferson Lab used polarized sohd NH3 and ND3 
targets ^^'^^ based on the same technique of dynamic nuclear polarization. We refer 
interested readers to Ref. for details on these targets. 

4.1.3. Recoil polarimeters 

Focal plane polarimeters were standard equipment at intermediate energy proton 
accelerators, such as LAMPF ^o, TRIUMF ^i, SATURNE ^2, and PSI ^3. Sensi- 
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tivities of spin observables to small amplitudes were demonstrated by experiments 
carried out at these facilities using focal plane polarimeters. In recent years, proton 
recoil polarimeters have been employed in a number of experiments at electron ac- 
celerator facilities, such as the MIT-Bates Laboratory the Mainz laboratory 
and the Jefferson Laboratory 56.57,58 

Proton polarimeters are based on nuclear scattering from an analyzer mate- 
rial like carbon. The proton-nucleus spin-orbit interaction results in an azimuthal 
asymmetry in the scattering distribution which can be analyzed to obtain the pro- 
ton polarization and spin orientation. A typical proton recoil polarimctcr consists 
of two front detectors to track incident protons, followed by a carbon analyzer and 
two rear detectors to track scattered particles. Recoil neutron polarimeters which 
have been used in experiments to extract neutron electric form factor are based on 
the same physics principle and are constructed in similar ways. 
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Fig. 6. World data on the proton electric and magnetic form factor ratio, -qT' > "■^ function of 

M 

. The high precision data from polarization transfer measurements at Jefferson Lab by Gayou 
et al. ®* and by Jones et al. ^® plotted with statistical uncertainties only, show different Q^- 
dependence of the form factor ratio from those obtained by the Rosenbluth method from unpolarized 
cross-section measurements, at large values of . 



4.2. Polarized Electron- Proton Elastic Scattering 

While precise information on and is important for understanding the un- 
derlying electromagnetic structure of the proton, it is also very interesting to study 

the ratio of these two form factors, as a fimction of Q^. Any dependence in 

the proton form factor ratio would suggest different charge and magnetization spa- 
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tial distributions inside the proton. Double polarization experiments allow precise 
measurement of the proton form factor ratio directly. Furthermore, by combining 
polarization measurements with differential cross-section measurements, one can 
determine and separately with significantly improved precision compared 
to the conventional Roscnbluth separation technique. 

In the onc-photon-exchange Born approximation, the elastic scattering of lon- 
gitudinally polarized electrons from unpolarized protons results in a transfer of 
polarization to the recoil proton with only two nonzero components: Pj perpendic- 
ular to, and Pi parallel to the proton momentum in the scattering plane . The 
form factor ratio can be determined from a simultaneous measurement of the two 
recoil polarization components in the scattering plane: 

PtE + E' 9 

where E and E' are the incident and scattered electron energy, respectively, and 
6 is the electron scattering angle. 

Although no dramatic Q^-depcndcncc in the proton form factor ratio has been 
observed from unpolarized measurements, new data from a polarization transfer 
experiment which measured this ratio directly shows very intriguing behavior at 
higher Q^, i.e., G^ falls off much faster than G^ as a function of Q^. The intriguing 

result on at high values from Jefferson Lab has motivated much interest 

on this subject both experimentally and theoretically. A more recent polarization 
transfer experiment at Jefferson Lab showed that the downward trend continues 
out to Q"^ = 5.6 (GeV/c)^. Prior to the Jefferson Lab experiments, the recoil 
polarization technique was employed at the MIT-Bates laboratory and it has 
also been used more recently at Mainz Fig. 6 shows the world data on the proton 
form factor ratio with statistical uncertainties only from these recoil polarization 
measurements, as well as those obtained from unpolarized experiments. The new 
data from polarization experiments at Jefferson Lab show imprecedcnted precision 
compared to the data from unpolarized experiments. These new data demonstrated 
an intriguing dependence, dramatically different from the behavior observed in 
the unpolarized Roscnbluth separation experiments. Recently, a new experiment 
was completed at Jefferson Lab aiming at resolving such a discrepancy by using 
the so-called super Rosenbluth separation technique. Results from this experiment 
are anticipated in the near future. 

Alternatively, one can determine the proton form factor ratio by measuring the 
asymmetry of longitudinally polarized electrons scattered from a polarized proton 
target. The one-photon-exchange diagram for such spin-dependent elastic scattering 
is shown in Fig. 7. For longitudinally polarized electrons scattering from a polarized 
proton target, the differential cross section can be written as ^^: 

1 = ^ + ^^' 
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where E is the unpolarized differential cross section given by Eqn. 9, /i is the 
electron helicity, and A is the spin-dependent differential cross section given by: 

A = -(TMott/,7eLj2r«T' cos0*GS/ - 2^2t{1 + t)vtl' sm0* cos 0*G^,,GP ] , (20) 




Fig. 7. The one-photon-exchange diagram for spin- dependent electron-nucleon elastic scattering. 

where 6* and 0* are the polar and azimuthal proton spin angles defined with 
respect to the three-momentum transfer vector q and the scattering plane (Fig. 7), 
and vt' and vtl' are kinematic factors 

The spin-dependent asymmetry A is defined as: 



A 



h- ■ 



(21) 



where a'^^ denotes the differential cross section for the two different helicities 
of the polarized electron beam. The spin-dependent asymmetry A can be written 
in terms of the polarized and unpolarized differential cross-sections as: 

^ — / '/Tl I -t- T 17Irnr / Sm H ' COS It) I t*. ,( t'„ 

(22) 



^_A_ 2TVT'Cose*Glj 

^ S ^ 



2y/2HTT7jvTL' sine* cos 0*G5^G^ 



{1 + t)vlG''/ + 2tvtGI/ 

The experimental asymmetry A^^p, is related to the spin-dependent asymmetry of 
Eqn. 22 by the relation 



Aexp — PbPfA , 



(23) 



where Pb and Pt are the beam and target polarizations, respectively. A deter- 
mination of the ratio 7# , independent of the knowledge of the beam and target 
polarization, can be precisely obtained by forming the so-called super ratio 



_ 2^2t(1 + t)vtl' sin0l cos (/.^G^.^G^ 

^2 ~ 2tvt' cos 6I^G^/ - 2y'2T{l + T)vTL' sin6i* cos ^i^G^^G^ 



Ai 2tvt' cos eiG\/ 



(24) 
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where Ai and A2 are elastic electron-proton scattering asymmetries measured at the 
same value, but two different proton spin orientations corresponding to {Oi,(j)l) 
and (^2!'/'2)) respectively. For a detector configuration that is symmetric with 
respect to the incident electron momentum direction, Ai and A2 can be measured 
simultaneously by forming two independent asymmetries with respect to cither the 
electron beam helicity or the target spin orientation in the beam-left and beam- 
right sector of the detector system, respectively. Thus, the proton form factor ratio 
can be determined with high precision using this technique, and an experiment is 
currently planned at the MIT-Bates Laboratory Such an experiment is essential 
because it employs a different experimental technique that has different systematic 
uncertainties than those from recoil proton polarization experiments, and those 
from Rosenbluth and super Rosenbluth measurements. 

4.3. Polarized Quasielastic Electron-Nucleus Scattering and Neutron 
Electromagnetic Form Factors 



4.3.1. Experiments with Deuterium, Targets 

Measurements of the neutron electric form factor are extremely challenging be- 
cause of the lack of free neutron targets, the smallness of the G^, and the dom- 
inance of the magnetic contribution to the unpolarized differential cross-section. 
A promising approach to measure is by using polarization degrees of freedom. 
For coincidence elastic scattering of longitudinally polarized electrons from "free" 
neutrons, n{e,e'n) process, the recoil neutron polarization is given by 



hPe E + E' 



eV^(TT7)G^^an2(ie) (25) 



P. = 



hPe 



.(_2eV7(rT7)Gl;G^)tan(i0), (26) 



where P^ and P^ are the recoil neutron polarization component along the recoil 
neutron momentum direction, and transverse to it, respectively, in the scattering 
plane. Here Pg denotes the absolute value of the electron polarization and h = ±1 
denotes the electron helicity. In the case of a "free" neutron target, the normal 
component Py is zero due to time and parity invariancc. However, Py is typically 
not zero because an effective neutron (nuclear) target must be employed and final 
state interactions play a role. Thus, by forming the ratio of the unpolarized 
cross-section as well as the electron polarization cancels out, providing sensitive 
experimental access to the neutron form factor ratio. 



Pr -2M 



Pz E + E' 



tan(-( 



-1 



Eqn. 27 is only valid in the PWIA picture and corrections need to be applied 
when effective neutron targets are used. The deuteron is the simplest effective 
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neutron target and it is also the most effective neutron target for recoil neutron po- 
larization measurements. Correcting for nuclear binding effect, relativity, reaction 
mechanisms such as FSI and MEC effects is essential in order to extract reliable 
information on from quasielastic d{e,e'n) reaction. Furthermore, precise in- 
formation on is crucial in order to extract because it is the ratio to 
first order that is determined from such a recoil neutron polarization measurement. 
Experiments with longitudinally polarized electron beams and recoil neutron po- 
larimeters have been car ricd out at MIT-Bates and Mainz *53.64 ^j^g relatively 
low region, and has been extracted from the c?(e, e'n) process, using the 
state-of-the-art two-body calculations by Arenhovel Most recently, such an ap- 
proach has been employed at Jefferson Lab up to a value of 1.5 (GeV/c)^ 

Alternatively, one can employ a polarized deuteron target to probe the neu- 
tron electric form factor by the d{e, en) reaction. The scattering cross-section for 
longitudinally polarized electrons from a polarized deuteron target for the d{e, en) 
reaction can be written as 

S = So{l + PfA^ + P2AI + h{A, + PtAX^ + P^AX^)] , (28) 

where Sq is the unpolarized differential cross section, h the polarization of the 
electrons, and Pf (Pf) the vector (tensor) polarization of the deuteron. Ae is the 
beam analyzing power, A^^'^ the vector and tensor analyzing powers, and A^j^^ 
the vector and tensor spin-correlation parameters. The polarization direction of the 
deuteron is defined with respect to the three-momentum transfer vector, q. The 
vector spin-correlation parameter A^^ contains a term representing the interference 
between the small neutron electric form factor and the dominant neutron magnetic 
form factor, when the target spin is perpendicular to the q vector direction. Thus, 
the spin-dependent asymmetry (defined with respect to the electron beam helicity) 
from the d{e, en) reaction for vector polarized deuteron gives access to the quantity 

..... . -* 

7# to first order when the target spin direction is aligned perpendicular to q. Such 

experiments are extremely challenging since they involve both neutron detection 
and a vector polarized deuteron target. 

Recently, the neutron electric form factor was extracted for the first time 
from a d(e, e'n) measurement in which a vector polarized deuteron target from an 
atomic beam source was employed. Most recently, a d{e, e'n) experiment using a 
dynamically polarized solid deuterated ammonia target was carried out at Jefferson 
Lab and G^ was extracted at a value of 0.5 (GeV/c)^. The precision of G^ from 
these polarization experiments is significantly better than those by Platchkov et 
al. from the unpolarized electron-deuteron elastic scattering measurement. Thus, 
by using polarization degrees of freedom and coincidence detections of scattered 
electrons and recoil neutrons, the dominant neutron magnetic contribution and 
the proton contribution to the scattering process are suppressed and more precise 
information on the neutron electric form factor can be extracted. 

4.3.2. Quasielastic ^He{e,e') and ^He{e,e'n) Reactions 
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A polarized '^He nucleus is effective in studying the neutron electromagnetic 
form factors because of the unique spin structure of the "^He ground state, which is 
dominated by a spatially symmetric S wave in which the proton spins cancel and 
the spin of the ^He nucleus is carried by the unpaired neutron ^8,69^ Y\g, 8 shows 
the one-photon exchange diagram for longitudinally polarized electrons scattering 
off from a polarized ^He nuclear target at the quasielastic kinematics in the PWIA 
picture. For inclusive measurement, only the scattered electron is detected; for the 
coincidence ^He{e, e'n) measurement, both the scattered electron and the knockout 
neutron are detected. 




Fig. 8. The one-photon-exchange diagram for spin- dependent guasielastic electron-^ He scattering 
in the plane-wave impulse approximation picture. 



The spin-dependent contribution to the inclusive '^He(e, e') cross section is com- 
pletely contained in two spin-dependent nuclear response functions, a transverse 
response Rt' and a longitudinal-transverse response Rtl' '^^ ■ These appear in ad- 
dition to the spin- independent longitudinal and transverse responses Rl and Rt- 
These spin-dependent response functions Rt' and Rtl' can be isolated experimen- 
tally by forming the spin-dependent asymmetry A defined previously with respect 
to the electron beam helicity. In terms of the nuclear response functions, A can be 
written 

^ _ -{cosO*vt'Rt' + 2 sing* cos (I)*vtl'Rtl') ^29) 
vlRl + vtRt 

where the Vk are kinematic factors, and 9* and (jf are the target spin angles de- 
fined previously. The response functions Rk depend on and the electron energy 
transfer uj. By choosing 9* = 0, i.e. by orienting the target spin parallel to the mo- 
mentum transfer q, one selects the transverse asymmetry At' (proportional to Rt')', 
by orienting the target spin perpendicular to the momentum transfer q {9* — 90, 
4>* =0), one selects the transverse- longitudinal asymmetry Atl' (proportional to 
Rtl')- 
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Fig. 9. The neutron magnetic form factor data published since 1990, in units of the standard 
dipole form factor parameterization Gu, as a function of . The data by Gao et al. Xu et 
al. '^9.80 are from experiments employing polarized He targets discussed in the text, plotted with 
statistical uncertainties only. The rest are from experiments involving deuterium targets. The 
points of Anklin 94 and Gao 94 have been shifted slightly for clarity. 



Because the ^He nuclear spin is carried mainly by the neutron, Rt' at quasi- 
elastic kinematics contains a dominant neutron contribution and is essentially pro- 
portional to (G"|) , similar to clastic scattering from a free neutron. Unlike the 
free neutron, however, the unpolarized part of the cross section (the denominator 
in Eq. (29)) contains contributions from both the protons and the neutron in the 
^Hc nucleus at quasielastic kinematics. Therefore, At' is expected to first order to 
have the form (G^) V(a + 6(G^)^) in PWIA, where a is much larger than 6(G^)^ 
in the low region. While measurements of G^ using deuterium targets enhance 
the sensitivity to the neutron form factor by detecting the neutron in coincidence, 
a similar enhancement occurs in inclusive scattering from polarized ^He because 
of the cancellation of the proton spins in the ground state. This picture has been 
confirmed by several PWIA calculations ^-'^^-'^^^ more recent and more advanced 
calculation which fully includes FSI and most recently by Golak et al. in 
which both the FSI and MEG effects have been taken into account. 

The first experiments which investigated the feasibility of using polarized 
^He targets to study the neutron electromagnetic structure were performed at the 
MIT-Bates Linear Accelerator Center. Following these two experiments, the first 
measurement of G^ from •^He(e, e'n) was reported by Meyerhoff et al. and the 
neutron magnetic form factor at low was extracted for the first time ^® from 
inclusive measurement of the quasielastic transverse asymmetry, At' ■ Recently, a 
high precision measurement of At' was carried out in Hall A at Jefferson Lab at 
six quasielastic kinematic settings corresponding to central values of 0.1 to 0.6 



22 Nucleon Electromagnetic Form Factors 

(GeV/c)^. A state-of-the-art non-rclativistic Faddeev calculation was employed 
in the extraction of G%[ at the two lowest values of the experiment and a PWIA 
calculation was used in the extraction of G]^ for the remaining values of the 
experiment 

The Faddccv calculation treats the ''He target state; and the 3N scattciring states 
in the nuclear matrix element in a consistent way by solving the corresponding 3N 
Faddeev equations The MEG effects were calculated using the prescription of 
Riska , which includes tt- and p-like exchange terms. The Faddeev based formalism 
has been applied to other reaction channels and good agreement has been found 
with experimental results in particular with the recent NIKHEF data on Ay at 
Q2 _ Q iQ (GeV/c)^ from the quasielastic ^He(e, e'n) process and with the high 
precision '^He quasielastic asymmetry data in the breakup region at low from 
Jefferson Lab. This calculation, while very accurate at low Q^, is not expected to 
be sufficiently precise for a reliable extraction of from the "^He asymmetry data 
in a relatively higher region because of its fully non-relativistic nature. 

Fig. 9 shows data from the Jefferson Lab experiment '^^,80^ ^j^g earlier MIT-Bates 
experiment and data from deuterium experiments since 1990. The extracted 
values of at of 0.1 and 0.2 (GeV/c)-^ are in good agreement with previous 
measurements of Anklin et al. ^3,35 rpj^g Jefferson Lab data provide the first 
precision information on G^^ using a fundamentally different experimental approach 
than previous deuterium experiments. The extracted G^ values *° from the same 
experiment based on the PWIA calculation are also shown. While limitations 
exist in such an approach due to theoretical uncertainties, the results are in very 
good agreement with the recent deuterium ratio measurements from Mainz 35,36^ 
and in disagreement with results by Bruins et al. 

Enormous progress has been made in extracting G^ from inclusive quasielastic 
transverse asymmetry At' measurement in the low region due to the new preci- 
sion data and recent advances in three-body calculations. On the other hand, Atl' 
from quasielastic ■^He(e, e') process at low (Q^ < 0.3(Gey/c)^) is dominated 
by the proton contribution largely because of the smallness of G^ and the small 
non-S' wave part of the ^He ground state wave function. Thus, it is unreliable to 
extract information on G^ at low from the inclusive ^He{e, e') process. It 
is possible to go to higher Q'^ {Q^ > 0.3{GeV/c}^) to extract G^ with respectable 
accuracy from quasielastic ^He(e, e') mcasm-ement where the proton contribution 
to Atl' is under better control. It is more advantageous to determine the neutron 
electric form factor from the ■^He(e, e'n) reaction rather than the "^ife(e, e') reaction 
because the proton contribution to the asymmetry is minimized. 

The experimental asymmetry for the coincident reaction ^He{e, e'n) at the 
quasielastic kinematics can be expressed as follows in PWIA: 



2y^HjTTjta.n{e/2)G%Glj sm{0*) cos(0*) 
Gl^ + Gl/ij + 2r(l + r) tan2(0/2)) 
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p p ^ 2rVl + r + (1 + r)2 tan^(g/2) t.an(g/2)GS/ cos(r ) 
' " Gl;2 + G^2(^^2T(l + T)tan2(e/2)) ^ ' 

Here Pe is the electron polarization. P„ is the neutron polarization, D is an overall 
dilution factor which contains dilution from (possible) unpolarized neutrons in the 
target and dilution from background neutrons generated in {p, n) reactions, e.g. in 
shielding walls. Eqn. 30 shows the obvious sensitivity to in the longitudinal- 
transverse interference term. Therefore, by aligning the target spin perpendicular 
to q {A = Aperp), i.e. choosing 0* = 90°, and </>* = 0° the above equation can be 
rewritten in the following form {G^^ w 0): 

^„ ^ Aperp G^(r + 2T(l + r)tan^(^/2)) 
^ PePnD 2v'r(l+T)tan(0/2) ^ ' 



Aligning the target spin parallel to q (A = Ai,,,,,,) reduces Eqn. 30 to (G^;^ » 0): 



Along— -Tei-n-L' i , r>/i , _\ ^ ., , 9 /n /o\ " ^'^^1 



^2^1 + r + (1 + r)2 tan^(6'/2) ta n(6l/2) 
1 + 2(1 + t) tan2(0/2) 

This equation is completely independent of the neutron form factors and serves 
as an excellent calibration reaction. Thus, one can combine the above two equations 
and obtain 

-A 



GI = Jt + r(l + r) tan2(0/2)-i^G^. (33) 

■^long 

Precise information on GJ^ is a priori requirement before one can extract G^ 
from such double polarization measurements. The PWIA picture described above 
is obviously over-simplified. To extract reliably, corrections need to be ap- 
plied, in particular the correction for the FSI effect, which has been carried out 
recently using the Faddeev approach , There is also technical difficulty in- 
volved with detecting knocked out neutrons. Such experiments have been carried 
out at Mainz and NIKHEF in recent years. The first measurement on from 
^Hc(e, e'n) was reported by Meycrhoff et al. in which a high pressure polarized 
^He target achieved by the metastability-exchange optical pumping technique and 
the compression method was employed. More recent Mainz measurement by Becker 
et al. and Rohe et al. using the same technique show much improved statistical 
accuracy. The FSI correction has been applied in the extracted G^ value from 
the measurement by Becker et al. which is shown in Fig. 10. While the G^ 
shown in Fig. 10 by Meyerhoff et al. and Rohe et al. ^® do not include corrections 
for FSI, the FSI effect is expected to be small around a value of 0.6 (GeV/c)^ 
and the FSI effect is expected to increase the value of G^ by Meyerhoff et al. to a 
value of around 0.05. Also shown in Fig. 10 are the recent published G^ data from 
double-polarization deuterium experiments discussed previously, the extracted G^ 
values from the deutcron quadrupole form factor data by Schivilla and Sick and 
the Galster parameterization 
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Also shown are the extracted 



values from the deuteron quadrupole form factor data by Schivilla and Sick . 
parameterization is also shown. 
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Enormous progress has been made in the study of the neutron electromagnetic 

form factor using polarization degrees of freedom, and in particular with spin- 
polarized ■^He nuclear targets. However, it is still a very challenging task for both 
experimenters and theorists. The non-relativistic Faddeev approach has been well 
tested in the low momentum transfer and low energy transfer regime and reliable 
information on the neutron electric and magnetic form factor can be extracted in 
the low region. The development of fully relativistic three-body calculations 
is required to extract the neutron electromagnetic form factors reliably in the rel- 
atively higher region, which is an extremely difficult task. Recently, Kamada 
and collaborators solved the relativistic Faddeev equations for the three-nucleon 
bound state with a Lorentz boosted two-nucleon potential which is introduced in the 
context of equal time relativistic quantum mechanics. Such a calculation opens the 
door to considering the relativistic Faddeev equations for three-nucleon scattering. 



5. Theoretical Progress 

While enormous experimental progress has been made on the subject of nucleon 
electromagnetic form factors, significant theoretical progress has also been made 
in recent years in understanding the nucleon electromagnetic structure from the 
underlying theory of QCD. The newly discovered Generalized Parton Distributions 
(GPDs) 92,93^ which can be accessed through deeply virtual Compton scattering 
and deeply virtual meson production, connect the nucleon form factors and the 
nucleon structure functions probed in the deep-in-elastic scattering experiments. 
The GPDs provide new insights into the structure of the nucleon, and provide 
possibly a complete map of the nucleon wave-function. 

QCD is the theory of strong interaction and has been extremely well tested in 
the high-energy region, i.e., in the perturbative QCD regime. Ideally, one should 
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calculate the nucleon electromagnetic form factors directly from QCD in the non- 
perturbative region to confront the data. Unfortunately, nobody knows how to solve 
QCD analytically in the non-perturbative regime. Lattice QCD calculations based 
on first principles of QCD, on the other hand, have shown much promise; and 
this field is evolving rapidly due both to improvements in computer architecture 
and to new algorithms. While pQCD gives prediction for the nucleon form factors 
in the perturbative region, QCD effective theories such as the chiral perturbation 
theory can in principle provide reliable prediction in the very low energy region. In 
between the low energy region and the pQCD regime, varioTis QCD-inspircd models 
and other phenomenology models exist. Thus, precision data in all experimentally 
accessible region is crucial in testing these predictions. The rest of the section is 
devoted to brief discussions of various theoretical approaches used to calculate the 
nucleon electromagnetic form factors. 

5.1. Scaling and Perturbative QCD 

QCD is an asymptotically free theory, i.e., the strong coupling constant as ^ 
as the inter-quark distance 0. Thus, one can solve QCD using the perturbation 
method in the limit of oo. Such an approach is the so-called perturbative 

QCD (pQCD), and specific dependence of the nucleon electromagnetic form 
factors can be obtained from the pQCD analysis. 

Brodsky and Farrar proposed the following scaling law for the proton Dirac 
(Fi) and Pauli form factor (^2) at large momentum transfers based on dimensional 
analysis: 

The helicity- flipping form factor, F2, is suppressed by ^ compared to the 
helicity-conserving form factor Fi. Such a prediction is a natural consequence of 
hadron helicity conservation. The hadron helicity conservation arises from the vec- 
tor coupling nature of the quark-gluon interaction, the quark helicity conservation 
at high energies, and the neglect of the non-zero quark orbital angular momentum 
state in the nucleon. In terms of Sach's form factors G^^ and G^, the scaling result 
predicts: constant at large Q . Such scaling results were confirmed in a 

short-distance pQCD analysis carried out by Brodsky and Lepage Considering 

the proton magnetic form factor at large in the Breit frame, the initial proton 
is moving in the z direction and is struck by a highly virtual photon carrying a 
large transverse momentum, q±'^ = Q^. The form factor corresponds to the ampli- 
tude that the composite proton absorbs the virtual photon and stays intact. Thus, 
the form factor becomes the product of the following three probability amplitudes: 
(i) the amplitude for finding the valence \qqq > state in the incoming proton; (ii) 
the amplitude for this quark state to scatter from the incoming photon producing 
the final three-quark state with coUinear momenta and (iii) the amplitude for the 
final three-quark state to reform a proton. The contribution of the proton Fock 
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oo because each additional constituent contributes a factor of ^ to the 



states other than the valence \qqq > quark state to the form factor is suppressed as 
amplitude in (ii) 

Based on this picture, Brodsky and Lepage obtained the following result within 
their short-distance pQCD analysis 
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Fig. 11. The scaled proton Dirac and Pauli form factor ratio: (upper panel) and p 

^1 ^1 
(lower panel) as a function of Q'^ in (GeV/c)"^. The data are from Jones et al. ®® and Gayou et 
al. shown with statistical uncertainties only. 
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where as{Q^) and A are the QCD strong coupUng constant and scale parameter, 
bnm and 7^,71 are QCD anomalous dimensions and constants, and e|| (-e||) is the 
mean total charge of quarks with helicity parallel (anti-parallel) to the nucleon's 
helicity. For protons and neutrons 

pP — 1 _pP — Op" — _p" — _1 

e|| - i, e|| - U,e|| - e,, - ^, 

based on the fully symmetric flavor-hclicity wave function. While the constants 
bnm and C are generally unknown for baryons, they are equal for protons and 
neutrons using isospin symmetry. For the proton electric form factor, one obtains 
similar result for the dependence in the 00 limit and as such the short- 

distance pQCD analysis predicts the same scaling law as the dimensional analysis 
for the proton form factors: ^fi- — > constant and %A — » constant. Recently, 
Belitsky, Ji and Yuan performed a perturbative QCD analysis of the nucleon's 
Pauli form factor i^2(Q^) in the asymptotically large limit. They find that the 
leading contribution to F2{Q'^) goes like which is consistent with the scaling 

result obtained by Brodsky and Farrar The recent Jefferson Lab data ^^'^^ 

QP 

on the proton electric and magnetic form factor ratio, up to a Q value of 
5.5 (GeV/c)^ show a strong dependence in the ratio: the electric form factor 
falls off much faster than the magnetic form factor. In the naive picture, such 

dependence may suggest that the proton charge is distributed over a larger 
spatial region than its counterpart, the magnetization. Fig. 11 shows data on the 
scaled proton Dirac and Pauli form factor ratio '^pp^ (upper panel) and (lower 
panel) from Jefferson Lab experiments ^^^^''^ as a hmction of together with various 
predictions. The dash-dotted curve is a new fit based on an improved vector meson 
dominance model (VMD) by Lomon ^"^ . The thin long dashed curve is a point- 
form spectator approximation (PFSA) prediction of the Goldstone boson exchange 
constituent quark model (CQM) The solid and the dotted curves are the CQM 
calculations by Cardarelli and Simula ''^'^ including SU(6) symmetry breaking with 
and without constituent quark form factors, respectively. The long dashed curve is 
a relativistic chiral soliton model calculation The dashed curve is a relativistic 
CQM calculation by Frank, Jennings, and Miller 

While the short-distance pQCD analysis predicts a constant behavior for the 

in the Q — > 00 limit, the data are in better agreement with the -pj?- scal- 
. 1 . . . .1. 

ing behavior. The data could imply that the asymptotic pQCD scaling region has 

not been reached or that hadron helicity is not conserved in the experimentally 
tested regime so far. Such an experimental observation is consistent with studies 
from polarized deuteron photodisintegration polarized neutral pion production 
^H{j,p)tt'^ ^^'^j and deuteron tensor polarization measurements ed — > ed How- 
ever, Brodsky, Hwang and Hill were able to fit the Jefferson Lab data using 
a form consistent with pQCD analysis and hadron helicity conservation by taking 
into account higher twist contributions. Belitsky, Ji and Yuan also suggest that 
one should take into account the leading non-vanishing contribution in looking at 
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the pQCD scaling behavior between Ff and F2 form factors. On the other hand, 
Ralston and Jain argue that the scaling behavior is expected from pQCD 
when one takes into account contributions to the proton quark wave function from 
states with non-zero orbital angular momentum. 

Following an earlier work by Frank, Jennings and Miller Miller ^^"^ recently 
used light front dynamics in modeling the nucleon as a relativistic system of three 
bound constituent quarks surrounded by a cloud of pions. While the pion cloud 
is important for understanding the nucleon structure at low momentum transfers, 
particularly in understanding the neutron electric form factor, quark effects domi- 
nate at large momentum transfers. Such a model is able to reproduce the observed 
constant behavior of as a function of and the ratio is predicted to be 
a constant up to a Q"^ value of 20 (GeV/c)^. 

5.2. Lattice QCD Calculations 




0.0 ^ ^ ^ 1 ' ' ^ ^ 1 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

(GeV^) 

Fig. 12. Fit to the lattice QCD data for the square of the proton charge radius as a function of 
pion mass squared based on chiral extrapolation . The extrapolated value at the physical pion 
mass is shown by the solid dot with the large error bar, and the star indicates the more favorably 
accepted experimental value ^. 

The non-perturbative nature of QCD at large distances, i.e. low momentum 
transfers, prevents the analytical approach in solving QCD. However, important 

conceptual and technical progress has been made over the years, especially in the 
last decade or so in solving QCD on the lattice. While full QCD calculations have 
been carried out in some cases, most of the lattice results obtained so far were carried 
out in the so-called quenched approximation in which the quark loop contributions, 
i.e. the sea quark contributions, are suppressed. Furthermore, at present and in 
the foreseeable future, the lattice calculations are only practical using rather large 
quark masses because of the limitations of the currently available computational 
power. As such, uncertainties in extrapolating lattice results to the physical quark 
mass are rather large, particularly with the naive linear extrapolation in quark 
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mass. Thus, the chahenge is to find an accurate and reliable way of extracting 
the lattice results to the physical quark mass. The extrapolation methods which 
incorporate the model independent constraints of chiral symmetry particularly 
the leading non-analytic (LNA) behavior of chiral perturbation theory and the 
heavy quark limit arc exciting recent developments. Recently, such an approach 
has been used * in extrapolating the lattice result on the proton root-mean-square 
(rms) charge radius to the physical pion mass, and at present the extrapolated 
value is in good agreement with the more favorably accepted experimental value ^. 
Currently, lattice calculations of tli(^ nucleon form factor with smaller quark masses 
are underway The anticipated new results together with the improved LNA 
extrapolation methods will provide more reliable lattice results on the nucleon form 
factors in the near future. 

5.3. Other Calculations of Nucleon Electromagnetic Form Factors 
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Fig. 13. Jefferson Lab data 56, S8 proton form factor ratio, as a function of Q 

together with calculations from various vector meson dominance models 112,113,116 



In this section, we review various models and calculations of the nucleon elec- 
tromagnetic form factors. While direct comparison between data on proton and 
neutron electromagnetic form factors and model predictions is important, we will 
restrict the discussion to those calculations which agree reasonably well with the 
trend of the recent Jefferson Lab data on the proton form factor ratio and then 
look closely to see how well they describe individual form factors. The reason is 
that the Jefferson Lab data can already eliminate some calculations based on the 
high behavior, thus limiting the number of candidate calculations. In this spirit, 
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we display the calculations and the Jefferson Lab data versus In(Q^) which, while 
showing the Jefferson Lab data and the agreement (or lack thereof) of the various 
models, accentuates the low region where lattice QCD and chiral perturbation 
theory are promising in making reliable predictions. 

One of the earlier attempts to describe the proton form factors was a scmi- 
phenomenological fit by lachello and co-workers based on a model in which 
the scattering amplitude is written as an intrinsic form factor of a bare nucleon 
multiplied by an amplitude derived from the interaction with the virtual photon 
via vector meson dominance (VMD). Various forms of the intrinsic bare nucleon 
form factor were used: dipole, monopole, eikonal. However, since this function 
is multiplicative, it cancels out in the ratio G^/Gm- '^^^ VMD amplitude was 
written in terms of parameters fit to data. Gari and Kriimpelmann extended 
the VMD model to include quark dynamics at large via perturbative QCD. 
Because of the freedom in fitting the parameters, this model was later able to 
describe a subsequent set of data from SLAC that reported the ratio fXpG^/G^ 
growing significantly larger than 1 at large values of Lomon ^"^ extended the 

Gari and Kriimpelmann model by including the width of the p meson and the 
addition of higher mass vector meson exchanges. Such a model has been extended 
further to include the a;'(1419) isoscalar vector meson pole in order to describe 
the new Jefferson Lab proton form factor ratio data ■''S.ss pjg shows the new 
proton form factor ratio data as a function of together with predictions from 
various VMD models discussed above. While these models have little absolute 
predictive power, once the high data have fixed the parameters, the approach 
to low is highly constrained. One can clearly observe that there is significant 
difference between these calculations in the low domain. 

Hohler et al. fit a dispersion ansatz to e — A'' scattering data. VMD con- 
tributions from p, ui, (f), p', and cj' were included and parametrized. Prom these 
fits to data available at the time, a ratio iipG^^/G\j is obtained which is in good 
agreement with the present Jefferson Lab data up to Q-^ 3 GeV-^, as shown in 
Fig. 14, at which point the fitted curve begins to rise in contradiction to the trend 
of the; most recent data from Jefferson Lab. In addition, the best fit results in 
a value for the rms proton radius of 0.84 fm, about 4% lower than the currently 
accepted value, indicating that the low Q"^ behavior is not quite correct. In recent 
years, these VMD/dispersion relation approaches have been extended to include 
chiral perturbation theory Ii8'ii9420,i2i,i22 Mergell et al. obtained a best fit 
that gave an rms proton radius near 0.85 fm, closer to the accepted value of 0.86 
fm, but could not do better when simultaneously fitting the neutron data. Hammer 
et al. included the available data on the form factors in the time- like region in 
the fit to determine the model parameters. These fits, also shown in Fig. 14. tend 
to under-predict the Jefferson Lab results at low and over-predict at higher Q^, 
not being able to account for the slope. The latest work by Kubis which was 
restricted to the low domain of - 0.4 (GeV/c)^, used the accepted proton rms 
radius of 0.86 fm as a constraint. The available data between Q-^ of 0.05 (GeV/c)^ 
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and 0.4 (GeV/c)^ do not provide a severe test of the model, since they are hmited 
by the uncertainties in the measurements of in this region. However, Kubis' 
results show a marked departure from the trend of the Jefferson lab data as 
increases, decreasing far too rapidly as shown in Fig. 14. This is not unexpected 
since their region of validity was claimed to be for < 0.4 (GeV/c)^. 
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117,118,119,122 



Several authors have calculated the proton electric and magnetic form factors us- 
ing various versions of the constituent quark model (CQM). Chung and Cocstcr 
Aznauryan and Schlumpf ^-^^ all used a relativistic CQM (RCQM) to calculate 
nucleon form factors in the range of - 6 (GeV/c)^. Chung and Coester 
and Aznauryan were able to reproduce the then available data on F^iQ"^) and 

(Q^) in the range up to about 2-4 GeV^, above which the model deviated 
from the data particularly for Aznauryan. Schlumpf produced good agreement 
with later SLAC data showing a rise in the ratio of /UpG^ /^^m which has now been 
shown to be contradictory with the recent Jefferson Lab results. More recent cal- 
culations have been made using the CQM in light front dynamics (LFCQM) by 
Cardarelli, Simula, Pace, and Salme 126,127,99 rpj^jg approach uses a one-body cur- 
rent operator with phenomenological form factors for the CQs and light-front wave 
functions which are eigenvectors of a mass operator which reproduces a large part 
of the hadron spectrum. The SU (6) symmetry breaking effects with and without 
the constituent quark form factor are also included. Their calculations reproduce 
the trend of the Jefferson lab results ^^'^^ reasonably well. 

Frank, Jennings and Miller have calculated the proton form factors in the 
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CQM (in light front dynamics using the model of Schlunipf ^^^) with the primary 
focus of investigating the medium modifications in real nuclei. Their results for the 
free proton are in reasonable agreement with the data to the highest studied 
so far at Jefferson Lab and predict a change in sign of at slightly higher Q"^. 

Recently, this calculation was extended by Miller ^'^'^ . 

Li^^* used a relativistic quark model (RQM) in which symmetry is required in 
the center-of-mass frame. This has the effect of adding additional terms to the 
baryon wave function. The original wave function obtained before the SU{6) sym- 
metry requirement is constructed from spinors of quarks with zero momentum. The 
additional terms generated by the SU{6) symmetry requirement are constructed 
from spinors of anti-quarks with zero momentum. Taken together with the original 
terms, these represent the inclusion of the sea quarks. The results of this calculation 
originally preceded the publication of the Jefferson Lab results and the model 
gives good agreement with the data. 

Ma, Qing and Schmidt calculated the nucleon electromagnetic form factors 
within a simple relativistic quark spectator-diquark model using the light-cone for- 
malism. Melosh rotations are applied to both quark and vector diquark. They 
results describe the available experimental data well as shown Fig. 15. Recently, 
Wagenbrunn, Boffi, Klink, Plessas and Radici calculated the neutron and proton 
electromagnetic form factors for the first time using the Goldstone-boson-cxchange 
constituent quark model. The calculations are performed in a covariant framework 
using the point-form approach to relativistic quantum mechanics, and are shown in 
Fig. 15. 

Another recent calculation of this type is from a perturbative chiral quark model 
(pChQM) by Lyubovitskij et al. . In this model, the effective Lagrangian de- 
scribes relativistic quarks moving in a self consistent field, the scalar component 
of which provides confinement while the time component of the vector potential is 
responsible for short range fluctuations of the gluon field configurations, leading to 
a cloud of Goldstone bosons (tt, iiT, ry), as required by chiral symmetry. The results 
of these calculations drop much more rapidly than the experimental data. They also 
drop off more rapidly than either the results from the cloudy bag model (CBM), 
which uses a similar model, or those from the RCQM of Simula and co-workers. 

Holzwarth has used a Skyrme/soliton model to calculate the proton form 
actors and included loop corrections His results agree very well with the Jef- 
ferson Lab measurements. This results from a significant deviation of below 
the "standard" dipole form. In fact, it is a general feature of these models that 
changes sign for Q'^ somewhere larger than 5 GeV^ (although exactly where is not 
well predicted), while G^ either does not change sign at all or does so at a very 
large value of . They also agree very well with measured G^ and G\i values in 
the low region and yield a proton rms radius of 0.869 fm, in close agreement 
with the currently accepted value. 
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In the MIT Bag Model, culminating in the so-called cloudy bag model (CBM) 
which includes the bare MIT bag plus a meson cloud, quarks are described as 
independent particles confined in a rigid spherical well. Such calculations have 
been performed by Lu and Thomas 132,133,134 rpj^g introduction of the pion cloud 
improves the static properties of the nucleon by restoring chiral symmetry and also 
provides a convenient connection to ttN and NN scattering. The result of these 
calculations is a ratio of G^^/G\f which decreases much more rapidly than indicated 
by the; Jefferson Lab data. However, it provides a considerable improvement, as 
expected, over the bag model without the pion cloud which predicts a drop to zero 
at — 1.5 GeV^ with a subsequent change of sign. The CBM calculation ^^"^ with 
a bag radius of 0.7 fm is also shown in Fig. 15, and the abrupt change of slope 
between 2 and 3 GeV^ is inconsistent with the newer JLab results by Gayou et 
al. 

In summary, the Soliton model calculation by Holzwarth , the extended VDM 
model by Lomon the relativistic CQM model calculation by Miller which 

is an improvement over its earlier version ^'^^ , and the relativistic quark spectator- 
diquark model calculation by Ma, Qing and Schmidt describe the new Jefferson 
Lab proton form factor ratio data very well. Fig. 16 shows the comparison between 
these models and the Jefferson Lab data ^^'^s. 

While it is important for a model to explain the observed dependence of 
the proton form factor ratio, more stringent tests of the model can be provided by 
comparing calculations with data on all four nucleon form factors over the entire ex- 
perimental accessible momentum transfer region. Fig. 17 shows the proton electric 
and magnetic form factor data published since 1970 obtained by the Rosenbluth 



34 Nucleon Electromagnetic Form Factors 



1.2 




Fig. 16. Jefferson Lab data ''^''^'^ on the proton form factor ratio, as a function of 

together with calculations by Lomon (dash- dotted), Miller ^"'^ (solid curve), Holzwarth 
(dotted curve), and Ma, Qing, and Schmidt (dashed curve). 
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Fig. 17. World data on proton electromagnetic form factor, and since 1970 in the unit 
of the standard dipole parameterization as a function of together with calculations by Lomon 
(dash-dotted). Miller (solid curve), Holzwarth ^^'^ (dotted curve), and Ma, Qing, and 
Schmidt (dashed curve). 
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method from unpolarized cross-section experiments, together with predictions from 
these models; Fig. 18 shows the corresponding neutron data and the model predic- 
tions. While most of the models are in good agreement with the proton magnetic 
form factor data, the departure of the model predictions from the data at high 
is not surprising. These models describe the trend of the Jefferson Lab proton 
form factor ratio data well, which are in disagreement with the proton form factor 
ratio formed by and G^ obtained by the Rosenbluth separation method from 
unpolarized cross-section measurements. As shown in Fig. 18, none of these mod- 
els are capable of describing the neutron electric and magnetic form factor data 
simultaneously in both the low and high region. While enormous theo- 
retical progress has been made in recent years, significantly improved theoretical 
understanding the nucleon electromagnetic form factor is urgently needed. 
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Fig. 18. World data on neutron electromagnetic form factor, and since 1990 as a function 
of together with calculations by Lomon (dash-dotted) , Miller ^^"^ (solid curve), Holzwarth 
(dotted curve), and Ma, Qing, and Schmidt ^-^^ (dashed curve). 



6. Future Outlook 

New, high precision data have been obtained in recent years from different elec- 
tron accelerator laboratories around the world, thanks to new high luminosity fa- 
cilities, and novel polarization experimental techniques. Future measurements ex- 
tending into even higher momentum transfer region will be carried out in the near 
future, particularly with the planned energy upgrade of CEBAF to 12 GeV at Jef- 
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ferson Lab . With this upgrade, the proton form factor ratio measurement can 
be extended to a value of 14 (GeV/c)^ Such an extension will provide more 
stringent tests of various pQCD approaches discussed above. Furthermore, the neu- 
tron electric form factor measurement can be extended to a value of 4 (GeV/c)^ 
using the ^He{e, e'n) reaction ^^"^ . and the neutron magnetic form factor measure- 
ment can be extended to a value of about 10 (GeV/c)^ employing the ratio 
technique of D(e,e'n)/D(e,e'p) This technique has been employed successfully 
recently using the CLAS detector at Jefferson Lab 

Precision measurements of nucleon form factor in the low momentum transfer 
region, especially in the limit of 0, which allows the determination of the 

nucleon electromagnetic radii, are equally important. This is a region where lat- 
tice QCD is anticipated to make reliable predictions in the near future and where 
effective field theories are expected to work. At MIT-Bates, the newly constructed 
BLAST detector will determine the nucleon electromagnetic form factor with high 
precision in the < 1.0{GeV/c)^ region using the novel experimental technique of 
scattering longitudinally polarized electrons in an electron storage ring from polar- 
ized internal gas targets with a large acceptance detector. The proton rms charge 
radius will be determined with much more higher precision from BLAST than from 
existing measurements from unpolarized electron-proton clastic scattering experi- 
ments. Such information will be sufficiently precise to allow high precision tests of 
QED from measurements of the hydrogen Lamb shift and equally importantly, will 
provide reliable tests of Lattice QCD calculations. 
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